Effects of Lovelock terms on the final fate of gravitational collapse: 
analysis in dimensionally continued gravity 
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We find an exact solution in dimensionally continued gravity in arbitrary dimensions which de- 
scribes the gravitational collapse of a null dust fluid. Considering the situation that a null dust 
fluid injects into the initially anti-de Sitter spacetime, we show that a naked singularity can be 
formed. In even dimensions, a massless ingoing null naked singularity emerges. In odd dimensions, 
meanwhile, a massive timelike naked singularity forms. These naked singularities can be globally 
naked if the ingoing null dust fluid is switched off at a finite time; the resulting spacetime is static 
and asymptotically anti-de Sitter spacetime. The curvature strength of the massive timelike naked 
singularity in odd dimensions is independent of the spacetime dimensions or the power of the mass 
function. This is a characteristic feature in Lovelock gravity. 

PACS numbers: 04.20.Dw, 04.20.Jb, 04.50.+h 



I. INTRODUCTION 

Black holes are one of the most fascinating objects in 
general relativity. It is considered that they are formed 
by the gravitational collapse of very massive stars or den- 
sity fluctuations in the very early universe. It has been 
shown that spacctimes necessarily have a singularity un- 
der physically reasonable conditions p]. Gravitational 
collapse is one of the presumable scenarios that singular- 
ities are formed. 

Since there is no way to predict information from sin- 
gularities, we cannot say anything about the causal fu- 
ture of singularities. In order for spacetimes not to be 
pathological, Penrose made a celebrated proposal, the 
cosmic censorship hypothesis (CCH) Q, that is, singu- 
larities which are formed in a physically reasonable grav- 
itational collapse should not be seen by distant observers, 
i.e., spacetimes are asymptotically predictable. (See Q.) 
This is the weak version of the CCH. The strong version 
of the CCH states that no observer can see the singu- 
larities formed by gravitational collapse 0- Although 
there is a long history of research on the final fate of the 
gravitational collapse 0, S 0, Q, IM S E3, E3, E3I , 

we are 

far from having achieved consensus on the validity of the 
CCH, which is one of the most important open problems 
in general relativity. (See |l3l | for a recent review.) 

Lately it has been of great importance to consider 
higher-dimensional spacctimes. Although there is no di- 
rect observational evidence of extra dimensions, funda- 
mental theories such as string/M-theory predict the ex- 
istence of extra dimensions. Over the past years, the 
braneworld model has also reinforced the study of higher- 
dimensional spacetimes [Tij . 

There exists a natural extension of general relativ- 
ity in higher dimensions, Lovelock gravity [T^ j. Love- 
lock Lagrangian is composed of dimensionally extended 
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Euler densities, which include the higher-order curva- 
ture invariants with a special combination so that the 
field equations are of second order . In four dimen- 
sions, such higher curvature terms do not contribute to 
the Einstein field equations since they culminate in total 
derivatives. The quadratic curvature terms are known 
as Gauss-Bonnet terms, which appear in the low energy 
limit of hcterotic string theory |17|. In Gauss-Bonnet 
gravity, static black hole solutions have been found and 
investigated in detail [Isl. Il9|| . 

Higher-order curvature terms come into effect where 
gravity becomes very strong. One of the present authors 
has discussed the gravitational collapse of a null dust 
fluid in Gauss-Bonnet gravity [2(J • It was found that the 
final fate of the gravitational collapse in five dimensions 
is significantly different from that in general relativity: a 
massive timelike naked singularity can be formed. Such 
a class of singularities has never appeared in the general 
relativistic case The purpose of the present paper is 
to analyze how higher-order Lovelock terms modify the 
final fate of gravitational collapse in comparison to the 
Gauss-Bonnet or general relativistic cases. 

A static vacuum solution in Lovelock gravity has been 
found 0,H3| and studied (2j|. In this case, the metric 
function is obtained from an algebraic equation, which 
cannot be solved explicitly in general. This is because, 
in D-dimcnsional spacetimes, Lovelock Lagrangian con- 
tains [(-D + l)/2] independent free parameters, where the 
symbol [x] is understood to take the integer part of x. 
This complexity makes the analysis in Lovelock grav- 
ity quite hard to tackle, and hence it is not obvious to 
extract physical information from the solution. Conse- 
quently, Bahados, Teitelboim and Zanclli have proposed 
a method that reduces the \{D + l)/2] free parameters 
to only two by taking a special relation between the co- 
efficients of Lovelock Lagrangian [24[ . An exact solution 
in the theory of gravity, called dimensionally continued 
gravity (DC gravity) . repres enting a static vacuum solu- 
tion, has been found [24], [251 • DC gravity is very effective 
in estimating the effects of higher-order Lovelock terms 
although it is a restricted class of Lovelock gravity. 
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In this paper, we extend our investigations on the grav- 
itational collapse of a null dust fluid in Gauss-Bonnet 
gravity [2lJ into DC gravity in order to see the effects 
of higher-order Lovelock terms on the final fate of grav- 
itational collapse. An exact solution describing gravita- 
tional collapse of a null dust fluid is given. We discuss the 
final fate of gravitational collapse of a null dust fluid in 
comparison with the results in Gauss-Bonnet gravity |20j| . 

The outline of this paper is as follows; in Sec. [HI DC 
gravity is summarized. An exact solution for a null dust 
fluid is presented in Sec. 1 1 1 1 1 Sec. IIVI is devoted to the 
discussion of naked singularity formations in the gravita- 
tional collapse of a null dust fluid. The strength of naked 
singularities is investigated in Sec. Our conclusions 
and discussion follow in Sec. IVII The discussions on the 
geodesies from singularities are given in Appendix^ In 
Appendix [Bj the scalar curvature quantities are calcu- 
lated for the estimation of singular nature. We discuss 
the general relativistic case in Appendix [U] so as to com- 
pare it with the DC gravity results. We use the unit 
conventions c = 1 throughout the paper. 



II. DIMENSIONALITY CONTINUED GRAVITY 

We first summarize DC gravity as a restricted class of 
Lovelock gravity. Then we give the field equations in DC 
gravity for odd and even dimensions. 

Lovelock gravity is a natural extension of general rel- 
ativity in D(> 3)-dimensional spacetimcs, whose action 
is given by |l5| 



/= C 



-D 



with 



- [D/2] 

/ Cd = K ^2 a pl] 
p=0 



(2.1) 



(2.2) 



I p = t ai ... aD 1l aia2 Ml* 2 "- 10 - 2 " 

J M 

Ae Q2p+1 A--- Ae a °, (2.3) 



where 1Z 



ah 



duj ab 



■ A Lu c b is the curvature two-form, 



e a is the orthonormal frame, ui is the spin connection 
which satisfies the torsion-free Cartan structure equa- 
tions and I m is the action of the matter field, respec- 
tively. In our notation, a p , b p , ■ ■ ■ , runs over all spacetime 
dimensions 0, • • • , D — 1. e ai ... aD is the Z?-dimensional 
Levi-Civita tensor in Minkowski space (eoi - (D-i) = !)■ 
The [(D + l)/2] real constants a p in the action Ij2.2|l have 
dimensions [length] "t^ -2 ^; otherwise they are arbitrary. 
A real constant k has units of the action and can be set to 
be positive without loss of generality, (/tai)" 1 is propor- 
tional to the D-dimcnsional gravitational constant with 
a positive coefficient, so that it is real and positive. 

The Lovelock action (|2.1|) is the special combination 
of curvature invariants in that the resulting field equa- 
tions do not include more than the second derivatives 
of the metric. They are composed of dimensionally ex- 
tended Euler densities [la. In D-dimcnsions, the first 



\{D+l)/2] terms contribute to the field equations; higher- 
order curvature terms will be the total derivatives. In 
even dimensions, the last term in the summation (j2.1(l . 
that is p = D/2, is the Euler density and consequently 
does not contribute to the field equations. The action 
(|2.1(l reduces to the usual Einstein-Hilbert action in four 
dimensions. 

Varying the action (|2.1|l with respect to the vielbein 
form e a , one obtains the field equations: 

[D/2] 

-*Y1 a ^ D ~ 2p)e ai ... OD ft 0102 A • • • 

p=0 

A n a2p - ia2p A e Q2p+1 A • • • A e" - 1 = Q aD , (2.4) 

where Q ao is the (D — l)-form representing the energy 
momentum tensor T^ u of the matter, 



1 



(D-iy. 



T ao bl e bl ... bD e b2 A 



Ae 



hi, 



(2.5) 



The spherically symmetric vacuum solution in Lovelock 
gravity was first found by Whitt [2lJ ■ This solution has 
been extended to spacetime as a product manifold M. w 
M 2 x IC D ~ 2 by Cai where KP~ 2 is the maximally 
symmetric spaces with constant curvature (D — 2)(D — 
3)k. Without loss of generality, k is normalized as +1 
(positive curvature), (zero curvature) and —1 (negative 
curvature), respectively. 

The vacuum static solution in Lovelock gravity in- 
cludes \{D + l)/2] arbitrary constants 0, [22], so that 
the solution are given by the roots of polynomial of 
\{D — l)/2] degrees, which cannot be written down in 
explicit form. In DC gravity, the arbitrary constants 
are reduced to two by embedding the Lorentz g roup 
SO(D — l, 1) into the larger AdS group SO(D- 1, 2) [Hf. 
The remaining two arbitrary constants are a gravita- 
tional constant and a cosmological constant. Accord- 
ingly, Lovelock gravity is separated into two distinct type 
of branches for odd and even dimensions. The special 
combinations of Lovelock coefficients are given by 



1 ( n-\ 



D — 2p \ P 



l~D+2p 



[~D+2p for D = 2n _ i 



for D = 2n, 



(2.6) 



where n = [(D + l)/2] and l/l 2 is proportional to the 
cosmological constant with a negative coefficient. From 
Eq. H2.6f) . we obtain 



(«ai) 1 = 



(D - 2)/[(n - 1)kI- d+ % for D = 2n - 1 



\/{nKl- D + 2 ), 



for D = 2n, 
(2.7) 



so that I must be real and positive both for D = 2n — 1 
and D = 2n in order for the Z?-dimensional gravitational 
constant to be real and positive. 
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In odd dimensions D = 2n — 1, Lagrangian £271-1 is 
given by |24j 

n-l 

C-2n-l =K^a p e ai ... a2n _ 1 ?J° ia2 
p=0 

A ... A 7^«2p-i02 P A e a2 p +1 A • • ■ A e 02 "- 1 (2.8) 
with l|2.6|l . We adopt units such that 

{D - 2y.n D - 2 Kr 1 = 1 forD = 2n-l, (2.9) 
where fi.D-2 is the area of (D— 2)-dimcnsional unit sphere 

2 7r (^- 1 )/ 2 



D-2 



r((L>-i)/2)- 



(2.10) 



Eq. (|2.8|l is a D-dimcnsional Chern-Simons Lagrangian 
for AdS group, whose exterior derivative is proportional 
to the Eulcr density in (D + l)-dimcnsions. Defining two- 
forms 



ab 



K" =K ab + T 2 e a Ae b , 



(2.11) 



one can find the equations of motion for odd dimensions 
D = 2n- 1, 

-I Ke ai ... a2n _Jl A- ■■ATI =Qa 2 „-i- 

(2.12) 

In even dimensions D = 2n, Lagrangian is given by 

„ 2 CL1CL2 £ K2ji- 1&2TI 

£ 2n = Ke ai ... a2n n A---ATZ , (2.13) 

where we adopt the units as 

2D(D - 2)!f] D _ 2K r 2 = 1 forD = 2n. (2.14) 

Due to the particular choice of the Lovelock coefficients 
(|2.(i|) . the Lagrangian 12.13( 1 is brought into the Born- 
Infeld type 



£■271 = KV \detTZ 



(2.15) 



The field equations for even dimensions D = 2n are then 
given by 



~ ai Q2 



-2nl Ke ai ...a 2n 7l A ■ ■■ 

A ^a 2 , l - 3 a2„-2 



Ae" 2 - 1 =Q a2 „- (2.16) 



III. NULL DUST SOLUTION 

Bahados, Teitelboim and Zanelli derived the static and 
spherically symmetric vacuum solution of DC gravity in 
D(> 3)-dimensions [24||. The extension to spacetime 
M « M 2 x JC D ~ 2 has also been considered O. The 
general metric on the spacetime M. can be written as 

= <H&g{g A B,r 2 ^ij), (3-1) 

where §ab is an arbitrary Lorentz metric on M 2 , r is a 
scalar function on M 2 with r — defining the bound- 
ary of M 2 and 7^ is the unit curvature metric on K, D ~ 2 . 



When r is not constant, it is shown in full Lovelock grav- 
ity that the vacuum solutions are described by the static 
solutions obtained by Cai [23|; i.e, the generalization of 
the Birkhoff 's theorem into Lovelock gravity holds [2(| . 

From this fact, when r is not constant, th e g eneral 
vacuum solution in DC gravity is obtained as H3 



ds 2 = 



-f(r)dt 2 +.r 1 (r)dr 2 + r 2 dT,% 



D-2 



(3.2) 



with 



f(r) 



k- (2A//r) 1 /(«-i) +r 2 r 2 , 
fc-M 1 /^-!) +r 2 r 2 , 



for D = 2n, 

for D = 2n - 1, 
(3.3) 



aS ^k,D-2- 



where we denote the line element of KP~ 

It is noted that the conventions of M are slightly dif- 
ferent from those in [3 H3 ■ The constant M has dimen- 
sions of length in even dimensions, while it is dimension- 
less in odd dimensions. This is owing to the definition of 
the present unit conventions. (See Eqs. 1(2.9(1 and 1(2.14(1 .) 
The solution l|3.2|l is asymptotically D-dimcnsional AdS 
spacetime because I 2 > 0. 

The causal structures in the spherically symmetric case 
(k = 1) are shown in [24j . The solution represents a 
black hole for M > 0, and represents the D-dimcnsional 
AdS spacetime for M = in even dimensions. In odd 
dimensions, meanwhile, the solution represents a black 
hole for M 1 ^™ -1 ' > fc, the D-dimcnsional AdS spacetime 
for M = and the naked singularity for M 1 /^™ -1 ) < k. 
In three dimensions, this is the so-called BTZ solution 
|27tl28| . in which there are no curvature singularities. In 
four dimensions with k = 1, Eq. I|3.2|) reduces to the 
Schwarzschild-AdS solution. 

We will see in five and six dimensions, the solution 1(3.3(1 
belongs to the special family of solution in Gauss-Bonnet 
gravity. In five and six dimensions, up to quadratic cur- 
vature terms have nontrivial effects on the field equations. 
The Lagrangian in Gauss-Bonnet gravity is 

L = -2A + R + a{R 2 - AR^R^ + R^ pa R^ pa ). 

(3.4) 

The vacuum static solution in Gauss-Bonnet gravity 01 
is given by the metric l|3.2|l with 



/ 



r 

25 



1 +4<S 



(3.5) 



where a = (D - 3)(D - 4)a, A = 2A/[(D - l)(D - 2)] 
and 771 is the integration constant. We can find easily the 
solution (|3.5|) reduces to the solution l|3.3|l with 1+45A = 
and 2a = I 2 in five and six dimensions. 

In four-dimensional spacetimes, we can obtain the 
Vaidya solution by replacing the mass parameter M in 
the Schwarzschild solution as an arbitrary mass function 
M(v), which specifies the flux of the ingoing radiation 
with the advanced time coordinate v. In the course of this 
heuristic procedure, we find an exact solution of DC grav- 
ity in D-dimcnsional spacetimes by taking M — > M(y) 
in the solution (|3.2|l . representing a radially ingoing null 
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dust fluid. The energy momentum tensor of a null dust 
fluid is written as 



T, 



(3.6) 



where p(v, r) represents the energy density of the null 
dust fluid and Z M is the null vector normalized as l^d a — 



-d r . Then we can find the metric of the null dust solution 



as 



ds 2 = -f(v, r)dv 2 + 2dvdr + r 2 dY? kD _ 2 (3.7) 



with 



f(v,r) 



k- (2M(«)/r) 1 /(«-i) +r 2 r 2 

k - M^VC"- 1 ) +r 2 r 2 , 

I 



for D = 2n, 
for D = 2n- 1. 



(3.8) 



We can confirm the solution l|3.7|l to satisfy the field equa- 
tions with the help of 



n 



n 



2M(v) 



l/(n-l) 



e"Ae r , (3.9) 



2{n — l) 2 r 2 y r 

1 ( 2M{v) ^ 
~2(n- l)r 2 V r 

M(v) (2M(v)\ X/(n - 1) eVAei 



v.: 



K 



2(n - l)rM(v) \ r 

1 f 2M{v) \ 1/{n - 1] 
2(n- l)r 2 V r J 

1 [2M{v)\ 1,{n - 1) , ■ 
e A e- 7 , 



-l+l/(n-l) 



for even dimensions and 



ri M(v)M(v) 



2(n - l)r 



n 



(3.10) 
e v Ae\ (3.11) 

(3.12) 



(3.13) 
(3.14) 



for odd dimensions, where the overdot denotes the differ- 
entiation with respect to the coordinate v, and i, j rep- 
resent the indices on KP~ 2 . The matter terms are given 
by 



(D-2) 



(3.15) 



From the field equations, the energy density of the null 
dust fluid is given by 



p(v,r) 



1 



■M 



(3.16) 



both in odd and even dimensions. M > is required 
in order for the energy density of the null dust fluid to 
be non-negative. The solution (|3.7() is a generalization of 
the Vaidya solution. Eq. I|3.7J) reduces to Eq. I|3.2J) with 
dv = dt + dr/ f(r) when the mass function M is constant. 



IV. NAKED SINGULARITY FORMATIONS 

In this and the next sections, we study the final fate 
of the gravitational collapse of a null dust fluid in DC 
gravity by use of the solution obtained in the previous 
section. In what follows, we restrict attention to the 
spherically symmetric case, i.e., k = 1, for comparison 
with the result in Gauss-Bonnet gravity [2£j]. Hereafter, 
we shall call the solution H3.7(l with k = 1 the DC- Vaidya 
solution. For the case of constant M, we shall call the 
solution the DC-BTZ solution. In the Vaidya solution, 
the mass function M(v) is the Misner-Sharp mass. In 
the generalized Vaidya solution in Gauss-Bonnet gravity, 
M{v) is not the Misner-Sharp mass but more preferable 
quasi-local mass [2(|. In this paper, we also adopt M{v) 
as the quasi-local mass in the DC- Vaidya solution. The 
odd-and even-dimensional cases are separately investi- 
gated below. 

Discussions with the fixed point method are also pre- 
sented in Appendix IA II 



A. Even dimensions 

We consider the situation in which a null dust fluid 
radially falls into the initial AdS spacctime (M(v) = 0) 
at v = 0. We set the mass function as a power law form 



M(v) = M v q 



(4.1) 



for simplicity, where Mq(> 0) and g(> 1) are constant. 
We can see easily from Eq. (I3.16|) that a central singu- 
larity appears at r = for v > 0. Let us discuss more 
specifically the singular nature of v = r = 0. 

The future-directed outgoing radial null geodesies obey 



— = -f 

dv ~ 2 ; ' 



(4.2) 



so a trapped region exists in / < 0, which is given by 

2M{v) = 2M v q > r(l + TV)™" 1 , (4.3) 

where the equality sign holds at the marginally trapped 
surface / = 0. Thus the singularity at r = 0, v > is in 
the trapped region and only the point v = r = has the 
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possibility of being naked. Along the trapping horizon, 
that is, the trajectory of the marginally trapped surface, 
we have 



ds 2 = 



AqM v 



9-1 



(1 + /- 2 r 2 )"- 2 (l + (2n - l)l- 2 r 2 ) 



dv 2 , (4.4) 



so that it is spacelike for v > 0, r = 0. 

We try to find the radial null geodesies emanating from 
the singularity in order to determine whether or not the 
singularity is naked. 

It is shown that if a future-directed radial null geodesic 
does not emanate from the singularity, nor does a future- 
directed causal (excluding radial null) geodesic. (See Ap- 
pendix for the proof.) So we consider here only the 
future-directed outgoing radial null geodesies. 

Suppose the asymptotic form of the geodesies as v ~ 
Kir p around v = r = |3lj , where p and Ki are positive 
constants. Then the lowest order of Eq. I|4.2[l around 
v = r = yields p = 1, so that 



v ~ K\T. 



(4.5) 



Ki = 2 is obtained for q > 1, while Ki satisfies the 
relation 



(4.6) 



i- y[ i 



for q = 1. The condition for the existence of a positive 
Ki satisfying Eq. (|4.6() is 



n „ , 1 fn- 1 
< M < — 



An V n 



(4.7) 



and then K\ > 2 holds. 

Along the radial null geodesic arising from v = r = 
with asymptotic form (|4.5|l . the energy density for the 
null dust fluid (|3.16J) and the Kretschmann scalar T\ = 
R^vpaR^ 9 " diverge for r — ► as (see Appendix lB| 



p = 0(l/r 



D-q-l 



and 



Ii = 0(1/ 



,4(l-(9-l)/(£>-2)) 



(4.8) 



(4.9) 



respectively, so they are singular null geodesies for 1 < 
q < D — 1. Thus, the spacetime represents the forma- 
tion of a naked singularity for 1 < q < D — 1 with any 
Mo(> 0) and for q = 1 with Mq satisfying the condition 
(|4.7|l . In the limit of n — > oo, the right hand side of Eq. 
(|4.7|) goes to zero; i.e., the formation of naked singularity 
is less plausible as the spacetime dimensions are higher. 
Hereafter we consider the case where 1 < q < D — 1, in 
which a naked singularity appears at v = r = 0. 

The central singularity v — r = is then at least 
locally naked. Now we consider the structure of the naked 
singularity. Expanding Eq. (|4.3[) around r = 0, we have 



(2M ) 



(4.10) 



For q > 1, there exists a spacetime region W which is 
both the past of the trapping horizon and the future of 



the future-directed outgoing radial null geodesic 7 which 
behaves as Eq. (|4.5|) near v = r = 0. Such a region also 
exists for q = 1 because 1 > 2MqKi holds from Eq. (|4. 
and Ki > 2. Because the trapping horizon is spacelike 
for v > and r > 0, the past-directed ingoing radial null 
geodesic £ emanating from an event in U never crosses 
the trapping horizon. Also £ never crosses 7 except for 
v = r = because they are both future-directed out- 
going radial null geodesies. Consequently, C inevitably 
reaches the singularity at v = r = 0. Since U is an 
open set, we then conclude that there exist an infinite 
number of future-directed outgoing null geodesies ema- 
nating from the singularity at v = r = 0. Such geodesies 
should correspond to the solution of Eq. i|A3(l with 77 = 
or i] = 00 in Appendix IA II On the other hand, the 
future-directed ingoing radial null geodesic reaching the 
singularity v = r = is only v — 0. Therefore, it is con- 
cluded that the singularity v = r = has an ingoing-null 
structure. 

When M fails to satisfy the condition 1)4. 7fl for q = 1 
in D > 4 (n > 2) dimensions, more detailed analyses are 
needed in order to determine whether the final fate of 
gravitational collapse is a black hole or a naked singu- 
larity since there might exist null geodesies which do not 
obey the power-law form asymptotically. 

In order to sec whether the singularity is globally 
naked, we consider a very simple situation in which the 
null dust fluid is switched off at v = Vf > 0. The so- 
lution is described by the DC-BTZ solution for v > 17, 
which is static and asymptotically AdS spacetime. The 
DC-BTZ solution is joined with the /^-dimensional AdS 
spacetime for v < by way of the DC-Vaidya solution 
for < v < Vf. The situation is depicted in Fig. ^ 

If a null geodesic emanating from the singularity v = 
r = reaches the surface v = Vf in the untrapped region, 
it can escape to infinity; then, the singularity is globally 
naked. If we take Vf to be sufficiently small, the singu- 
larity can be globally naked. The singularity in the outer 
DC-BTZ spacetime is spacelike in even dimensions [24^ . 
Thus, the Penrose diagram of the gravitational collapse is 
drawn in Fig. [21 for the globally naked singularity forma- 
tion. Of course, the locally naked singularity formation 
is also obtainable. 



B. Odd dimensions 

As in the even-dimensional case, we consider the sit- 
uation in which a null dust fluid radially falls into the 
initial AdS spacetime, and we choose the mass function 
as the power- law form M(v) = Mqv q . Then we can see 
from Eq. I|3.16fl that a singularity develops at r = for 
v > 0. As shown below, the point v — r = may also be 
singular. It is noted that the singularity with r = and 
v > is also interpreted as a conical sing ularity; that is, 
a deficit angle exists at r = 0. (See |23| for the conical 
singularity in the BTZ solution.) 

A trapped region / < is given by 



M(v) = M v q > (1 + rv) 



(4.11) 



where the equality holds at the marginally trapped sur- 
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cut horizon 




FIG. 1: Illustration of the gravitational collapse of a null dust 
fluid in even dimensions. A central singularity forms at v = 
r = 0. The DC-BTZ solution is attached to D- dimensional 
AdS spacetime via the DC-Vaidya solution. Zigzag lines de- 
note the central singularity. 




FIG. 2: The figure shows a Penrose diagram for a globally 
naked singularity formation. Infinity consists of a timelike 
surface J?. EH and CH denote the event horizon and the 
Cauchy horizon, respectively. 



face / = 0. Along the trapping horizon, we have 



ds 2 



M ql 2 v 



r(n - 1)(1 + l- 2 r 2 ) n ~ 2 



dv 2 



(4.12) 



so that it is spacelike for r > 0, v > 0. From Eq. i|4.11[l . 
the singularity r = might be naked for < v < vah, 
where vah = . Following the argument in Appendix 
IA 21 we consider only the future-directed outgoing radial 
null geodesies from the singularity. Suppose the asymp- 
totic form 



v ~ vq + K 2 r p 



(4.13) 



near r = with < vq < vah, where p and are 
positive constants. From the lowest order of Eq. 1)4.20 . 
we obtain 



1 



(4.14) 



for < vq < vah, while there are no null geodesies with 
the asymptotic form Eq. (|4.13|) for v = vah- 

Along the radial null geodesies from the singularity 
v = r = asymptoting to Eq. i|4.14[l . the Kretschmann 
scalar diverges as 



1 1 =0(1/^(1-9/(13-1))) 



(4.15) 



for D > 5. In the meanwhile, it is finite for D = 3. 
In fact, the scalar quantities for the DC-Vaidya solu- 
tion with D — 3 arc the same as those for the three- 
dimensional AdS spacetime, which are constant every- 
where; X\ = I2 = 12l~ 4 , R = —6l~ 2 . (The curvature in- 
variants X\ , X2 and R for the metric l|3.7|l are calculated 
in Appendix FBI) However, the energy density of the null 
dust fluid H3.16|) diverges along the null geodesies (|4.14|) 
as 



0(l/r 



,D-q-l 



(4.16) 



for D > 3. As a result, they are singular null geodesies 
for 1 < q < D - 1. 

On the other hand, the radial null geodesies from the 
singularity r — 0, < v < vah asymptoting to Eq. 
(|4.14(l arc singular null geodesies for any q{> 1). Along 
them, the Kretschmann scalar and the energy density 
of the null dust fluid diverge for D > 5 and D > 3, 
respectively. 

We can consider the situation as in the even- 
dimensional case; that is, the null dust fluid is turned 
off at some finite time v = Vf > in order to see whether 
the singularity is globally naked. Fig. [21 illustrates the 
situation in which the DC-BTZ solution for v > Vf is con- 
nected with the D-dimensional AdS spacetime for v < 
via the DC-Vaidya solution for < v < Vf . 

The structures of the singularity in the outer DC-BTZ 
spacetime in odd dimensions are spacelike, null and time- 
like for Vf > uah, Vt = vah and Vf < vah, respectively 
[24^ . Now we consider the structure of the naked singu- 
larity in the null dust region. The naked singularity with 
< v < vah is both an endpoint of the future-directed 
ingoing radial null geodesic and the initial point of the 
future-directed outgoing radial null geodesic, therefore 
we conclude that it is timelikc. 

For < Vf < vah, the singularity is always globally 
naked. In the case of Vf > vah, if we take vt as suffi- 
ciently close to vah , a null ray emanating from the singu- 
larity reaches the surface v = Vf in the untrapped region, 
so that it can escape to infinity; then, the singularity is 
globally naked. As in the even-dimensional case, globally 
naked singularities can be formed in odd dimensions. 

In the case of < vt < vah, the Penrose diagram of 
the gravitational collapse is shown in Fig. 0] The possible 
Penrose diagrams are depicted in Fig. [S]for Vf = vah 
and in Fig. Elfor Vf > vah- From the present analysis, 
however, it is not clear whether there exists a null portion 
of the naked singularity at r = for v = vah- 
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V = » = 



FIG. 3: Gravitationl collapse in odd dimensions for Vf > vah 
(left) and < Vf < vah (right). The DC-Vaidya solution for 
< v < Vf is jointed with the DC-BTZ solution for v > Vf 
and the AdS spacetime for v < 0. 




FIG. 4: Penrose diagram in odd dimensions for < Vf < 
vah- The massive timelike naked singularity appears. The 
resultant spacetime is always globally naked. 

V. STRENGTH OF NAKED SINGULARITIES 

As have seen in the preceding section, naked singular- 
ities can be formed in DC gravity. In this section, wc 
investigate the strength of naked singularities along the 
radial null geodesies with the asymptotic form l|4.5|l in 
even dimensions and (|4.14() in odd dimensions. We de- 
fine 

V' = R^k v , (5.1) 

where fc M = dx^/dX is the affincly parametrized tangent 
vector of the future-directed outgoing radial null geodesic 
with an affine parameter A. Wc evaluate the strength of 
naked singularities by the divergent behavior of ip in the 
neighborhood of the singularities. 

We consider the radial null geodesic with the tangent 
vector fc M , which emanates from a singularity at A = 
0. In four dimensions, the strong curvature condition 




(a) (b) 



FIG. 5: Possible Penrose diagrams in odd dimensions for vf = 
vah , which are globally naked. 




(a) (b) 



FIG. 6: Possible Penrose diagrams in odd dimensions with 
Vf > vah for the globally naked singularity formation. 



and the limiting focusing condition are satisfied along 
an afnnely parametrized geodesic if lim,\— >o X 2 ip > and 
limx^oXip > 0, respectively [33, H3 • However, these 
results have not been extended to the higher-dimensional 
case so far. 

Along the radial null geodesies emanating from the sin- 
gularity, we have 

k r = \fk\ (5.2) 
Then the straightforward calculation yields 

, = -M_2V,, (5 .3) 
where / is given by Eq. I|3.8fl . The geodesic equation for 
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the radial null geodesic is written as 



Even dimensions 



(5.4) 



As we will see below, the asymptotic form of the differ- 
ential equation (|5.4fl for the DC-Vaidya solution reduces 
to the following form: 



(5.5) 



where s(> 0) and a are constants. Then ip is expressed 
as 



il}~s(D-2)r a ~ 1 {k r ) 2 



(5.6) 



We suppose the asymptotic solution of Eq. H5.5J) around 
A = Aq as 



r(A)~£ 1 (A-A ) /3 , 



(5.7) 



where D\ and (3 are positive constants. The divergent 
behavior of ip is completely specified by the following 
two cases. 



Case 1 ) (3 ^ 1. Direct substitution of Eq. 
into Eq. (|5.5|) yields a = -1 and (3 = 1/(1 + s). 
Redefinition of the affine parameter enables us to set 
D\ = 1 and Ao = without loss of generality in such a 
way that A = corresponds to the central singularity. 
Then we have 



r(A)~A 1 /d+«). 



Therefore we obtain 



s{D-2) 2 
(1 + s) 2 ' 



and then it is concluded 



,. , 2 , s(D-2) n 
lim A 2 V = , ' > 0. 

A^O (1 + s) 2 



(5.8) 
(5.9) 

(5.10) 



Note that the right-hand side of Eq. (|5.10|) varies 
delicately with the scaling of the affine parameter. 

Case 2)0 = 1. We have to take into account the 
next leading term such as 



r(A) ~ D 2 {\ - Ao) + As(A - A ) 1+7 , 



(5.11) 



where 7 is a positive constant. The most dominant order 
of Eq. gives 7 = a + 1 and D 3 = -sD c * +2 /{(a + 

l)(a + 2)]. We set A = and D 2 = 1 by the redefini- 
tion of the affine parameter. Therefore the asymptotic 
solution is given by 

r(A) ~ A. (5.12) 

Hence we have 

^~s(L>-2)A Q -\ (5.13) 
it is then concluded 

lim A 1 ""^ = s(D - 2) > 0. (5.14) 

Discussions for the odd- and even-dimensional cases 
are made separately in the following subsection by the 
aid of the criterion obtained above. 



For q = 1, we find that the asymptotic form of Eq. 
()5.4(l near v = r = becomes 



with 



a\ r 



K Y -2 



(5.15) 



(5.16) 



where the constant K\ satisfies the equation l|4.6|l . This 
belongs to Case 1. We then obtain 

lim X 2 t/j = S '} D ~ 2 2 } > 0. (5.17) 

A^O (1 + Si) 2 V ^ 

We next examine the case where q > 1. In this case, 
the asymptotic form of the geodesic equation becomes 



dk r 



+ s 2 r (q ~ n)/{n ~ 1) {k r ) 2 =s (5.18) 



with 



n — 1 

This belongs to Case 2 with a = — 1 + (q — l)/(n — 1). 
Thus we obtain 

lim A 2(i-fe-i)/(B-2)ty = r D _ 2) > 0. (5.20) 

A^O 

The power of A depends on the power of the mass func- 
tion and the number of the spacetime dimensions. It 
takes values from for q = D — 1 to 2 for q = 1. The 
divergent behavior of ifj is summarized in Table [I] 



TABLE I: The divergent behavior of tp around v — r — 
along the radial null geodesies with the asymptotic form 14.511 
in even dimensions. 



q 


9 = 1 


1< q< D-l 


q = D- 1 


q> D- 1 





A" 2 


x -2(l-(q-l)/(D-2)) 


const. 






B. Odd dimensions 

The steps in the odd dimensions are quite similar to 
those in even dimensions. We will first consider the radial 
null geodesies from the singularity r = and < v < 
Vah with the asymptotic form (|4.14|l . The asymptotic 
form of the geodesic equation (|5.4|l is given by 



where the constant S3 is defined by 
2< Z (A/ ^) 1 /(«-i) 



0, 



S3 = 



(n- l)«o(l - (M i>g)V(«-i))2 



> 0. 



(5.21) 



(5.22) 
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This corresponds to Case 2 with a = 0. Consequently, 
we find 

lirn \ip = S3(D - 2) (5.23) 

for the singularity r = and < v < vah- It is em- 
phasized that the strength is independent both of the 
spacctime dimensions and the power of the mass func- 
tion. 

We finally consider the geodesies from the singular- 
ity v = r = 0. The asymptotic form of the radial null 
geodesies (|5.4I) becomes 

rlb r 

^L_ + S4r -l+«/(n-l)(jfer)2 s , 0j (5 24) 

with 

54 ^ (29M « )1/( "~%0. (5.25) 
n—1 

This corresponds to Case 2 with a = — 1 + q/(n — 1). 
Thus, we obtain 

lim A 2 * 1 -"/^- 1 ^' = sJD - 2) (5.26) 

for v = r = 0. The power of A takes the value from 
for q = D - 1 to 2(1 - 1/(D - 1)) for <? = 1. In 
particular, the strength of the singularity for r = and 
< v < vah is weaker than that for v = r = in 
the case of (D — l)/2 < q < D — 1, while in the case 
of q = (D — l)/2, their strengths are the same. The 
divergent behavior of ip around v = r = is summarized 
in Table HD 

TABLE II: The divergent behavior of ip at v — r = along the 
radial null geodesies with the asymptotic form 14.141 in odd 
dimensions. It is noted that, for the naked singularity with 
r — and < v < vah, if) oc A" 1 is satisfied independent of 
D and q. 



q 


9=1 


Kq<D-l 


q = D - 1 


q >D- 1 


4' 


^-2(1-1/(0-1)) 


A -2(l-q/(D-l)) 


const. 






VI. CONCLUSIONS AND DISCUSSIONS 

In this paper, wc analyzed the D(> 3)-dimcnsional 
gravitational collapse of a null dust fluid in DC grav- 
ity. We found an exact solution with the topology of 
M ~ M 2 x )C D ~ 2 , which describes the gravitational col- 
lapse. Applying this solution with spherical symmetry 
to the situation in which a null dust fluid radially in- 
jects into initially _D-dimensional AdS spacetime, we in- 
vestigated the effects of the Lovelock terms on the final 
fate of gravitational collapse. Our model is an exam- 
ple of dynamical and inhomogeneous collapse in Love- 
lock gravity and should provide better insight into the 
CCH in the context of higher-curvature gravity theo- 
ries. We supposed that (i) the power-law mass function 



M(v) = M v q , where M > and q > 1, and (ii) the null 
geodesies obey a power law near the singularity. 

We found that globally naked singularities can be 
formed in DC gravity. Furthermore, the final states of 
the gravitational collapse differ substantially depending 
on whether the spacetime dimensions are odd or even. 
This property is seen in the homogeneous collapse of a 
dust fluid in DC gravity [34l |3f|. It is also noted that 
odd-dimensional DC gravity is an exceptional case in the 
thin-shell collapse investigated in [3(|, only in which a 
naked singularity can emerge. 

A massless ingoing null naked singularity can appear 
in the even-dimensional case (D = 2n) for 1 < q < D — 1 . 
In the odd-dimensional case (D = 2n — 1), on the other 
hand, a massive timelike naked singularity is formed for 
any q(> 1). These naked singularities can be globally 
naked. In three dimensions, the singularity is not a cur- 
vature singularity, where the scalar curvatures are the 
same as those for the three-dimensional AdS spacctime 
although the energy density of a null dust fluid diverges. 
As a result, the formation of a naked singularity cannot 
be avoided in DC gravity, nor in general relativity. 

In DC gravity, massive timelike singularities appear in 
odd dimensions. When M(v) is intact in Eq. 1)4. 11|) . it is 
easily shown that this property is independent of the form 
of M(v) as long as M > and M (0) = 0. The formation 
of a massive timelike singularity in odd dimensions is con- 
sidered to be a characteristic feature in Lovelock gravity. 
Our results are consistent with those in general relativ- 
ity and in Gauss-Bonnet gravity. In general relativity, a 
massive timelike singularity appears in three dimensions, 
where the Einstcin-Hilbert term becomes first nontrivial. 
(See Appendix[ni) While in Gauss-Bonnet gravity, it ap- 
pears in five dimensions jwhere the Gauss-Bonnet term 
becomes first nontrivial |2Cj . 

In Sec. we investigated the strength of naked sin- 
gularities. We compare the strength by the divergent 
behavior of ip = R flv k fl k' / . The strength of the naked 
singularity at v = r = depends on D and q both in odd 
and even dimensions (see Tables HI and UTTl . In contrast, 
around the massive timelike singularity at r = and 
< v < vah in odd dimensions, ip diverges as A -1 , inde- 
pendent of D or q. Such a divergent behavior for the mas- 
sive timelike singularity can be seen as well both in gen- 
eral relativity in three dimensions (see Appendix Q and 
in Gauss-Bonnet gravity in five dimensions |2(| • Massive 
timclikc singularities diverging as A -1 might be salient 
features in Lovelock gravity. 

We have shown in this paper that naked singularity for- 
mation is generic for the power-law mass function in DC 
gravity. We can speculate from the results that naked sin- 
gularity formation will also be inevitable in full Lovelock 
gravity. It is also extrapolated that if massive timelike 
naked singularities appear in full Lovelock gravity, they 
would be developed only in odd dimensions and have a 
portion in which ip would diverge as A -1 regardless of the 
spacctime dimensions or the form of the mass function. 
Additional investigations are required in order to clarify 
these points. 



10 



Acknowledgments 

We would like to thank K. Maeda for continuous en- 
couragement. We would also like to thank T. Torii and 
U. Miyamoto for fruitful discussions and useful com- 
ments. This work was partially supported by the Waseda 
University Grant for The 21st Century COE Program 
(Holistic Research and Education Center for Physics Self- 
organization Systems) at Waseda University. 



APPENDIX A: GEODESICS FROM 
SINGULARITIES 

1. Existence and uniqueness 

We apply the fixed-point method for proving the exis- 
tence and uniqueness of future-directed outgoing radial 
null geodesies from singularities with asymptotic form 
P~5|) and 

We define new coordinates 9 and x as 



X=(v 



-v ) s 



(Al) 
(A2) 



where s is a positive constant and we only consider the 
region with v > vq. Then the geodesic equation (|4.2|) 
becomes 



dx sx 



with 



sx 



1 

2?7 



+ 



? 2 X 2/s 
2ry/ 2 



*(0,x) 



2r) 



(A3) 



(A4) 



where a new parameter rj (0 < r\ < oo) has been intro- 
duced. The form of <I> is given below for each case. 

If we choose the parameter rj and s to be t]q and So, 
respectively, with which $ is at least C 1 in x > 0, 6 > 0, 
then we can apply the contraction mapping principle to 
Eq. (|A3|) to find that there exists the solution satisfying 
9(0) = ?7o, and moreover that it is the unique solution of 
Eq. i|A3() which is continuous at x = 0. (See 0, [3(| for 
the proof.) 



a. even- dimensional DC-Vaidya solution 

For the even-dimensional DC-Vaidya solution, we have 
vq = and 



f2M v«- 1 \ 1/( - n - 1) 



V 9 

/2M oX (q ~ 1)/s 
{ 9 



l/(n-l) 



(A5) 



(A6) 



For q > 1, if we choose rj = ?/o = 1/2 and s = sq satisfying 



< s < 1, 
< s < 



(9-1) 



2(n- 1)' 



(A7) 
(A8) 



^ becomes at least C 1 in % > 0, 6 > 0, and consequently 
the existence of the solution with the asymptotic form 
(|4.5|l is shown. Unfortunately, this method cannot be 
applied to the case with q = 1. 



b. odd- dimensional DC-Vaidya solution 

For the odd-dimensional DC-Vaidya solution, we have 
$ = (M v q ) l ' {n - l \ (A9) 



[Mo(v +x 



lM<3ll/("-l) 



(A10) 



For < vq < vah j if we choose 77 = rjo = (1 — 
(M ^) 1/( "~ 1) )/2 and s = s satisfying < s < 1, * 
becomes at least C 1 inx>0,6>>0. For Vo = 0, on 
the other hand, if we choose r) = r)o = 1/2 and s = sq 
satisfying 



< s < 1 
< s < 



2(n- 1)' 



(All) 
(A12) 



^ becomes at least C 1 in x > 0, 9 > 0. Consequently, 
the existence of the solution with the asymptotic form 
(|4.14j) is shown. 



c. Vaidya-AdS solution 

We also study the general relativistic case. For the 
-D-dimensional Vaidya-AdS solution in general relativity 
(see Appendix EJ, we have 



$ = 



M v 



QD-3 ' 

M a (v + x 1/s ) q - D+3 

QD-3 



(A13) 
(A14) 



First we consider the case in D > 4 dimensions, in 
which we have vq = 0. For q > D — 3, if we choose 
7] = T)o = 1/2 and s = sq satisfying 



< s < 1, 



< s < 



?-(D-3) 



(A15) 
(A16) 



^> becomes at least C 1 in x > 0, 9 > 0, and consequently 
the existence of the solution with the asymptotic form 
(IC6|) is shown. Unfortunately, this method cannot be 
applied to the case with q = D — 3. 

Next we consider the three-dimensional case. For < 
v o < vah, if we choose r\ = 770 = (1 — MoVq)/2 and 
s = so satisfying < so < 1, \& becomes at least C 1 in 
X > 0, 9 > 0. For v = 0, on the other hand, if we choose 
rj = 770 = 1/2 and s = sq satisfying 



< s < 1, 
0<s < |, 

^> becomes at least C 1 inx>0,6'>0 



(A17) 
(A18) 
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2. On the null and causal geodesies 



and 



In Sec. IIVI we have analyzed only future-directed out- 
going radial null geodesies. The contraposition of the 
following theorem implies that it is sufficient to consider 
only the future-directed outgoing radial null geodesies 
in order to determine whether or not the singularity is 
naked. 

We will prove the following theorem: If a future- 
directed causal (excluding radial null) geodesic emanates 
from the central singularity, then a future- directed radial 
null geodesic emanates from the central singularity. The 
proof is similar to the four-dimensional case in |37j . 

We consider the DC-Vaidya solution, of which metric is 
given by (|3.7[) with k = 1. Without loss of generality we 
set the geodesies on the several equatorial plane thanks 
to the spherical symmetry. In the DC-Vaidya spacetime, 
the tangent to a causal geodesic satisfies 



dv, 



GG 



( dv\ dv dr L 2 
i \d\) + dXdX + ~r 2 = 



(A19) 



where A is an affine parameter, L 2 is the sum of the 
square of conserved angular momenta and e = 0,-1 for 
null and timelike geodesies, respectively. Then, at any 
point on such a geodesic, 



f ( —Y > 2— — 
J \d\) ~ dXdX 



(A20) 



with equality holding only for radial null geodesies. For 
the future-directed outgoing geodesies, this gives 



dv 2 



(A21) 



dr 



<^<0 

dr 



(A24) 



are satisfied in the untrapped and trapped regions, re- 
spectively, where the subscripts represent causal (ex- 
cluding radial null) geodesies and outgoing radial null 
geodesies, respectively. 

Let us consider the (r, w)-plane. The singularity is lo- 
cated at r — for v > 0. First we consider the sin- 
gularity with / < 0. By Eq. I|A24|) . the past-directed 
ingoing geodesies emanating from an event with r > in 
the trapped region cannot reach the singularity at r = 0. 
Therefore, there is no future-directed outgoing geodesic 
emanating from the singularity in the trapped region, 
i.e., the singularity with / < is censored. 

Next we focus on the singularity with < v < vah, 
where / > holds. Now suppose that v = vcG{ r ) ex- 
tends back to a central singularity located at (r, v) = 
(0,v s ), where v s satisfies < v s < Vah- There exists 
a portion of v = i | cg( j ") which is in the untrapped re- 
gion. Let p be any point on such a potion of v = vedf) 
and to the future of the singularity. Applying inequal- 
ity ljA23(l at p, we see that v = vrngM through p 
crosses v = ucgM from above and hence the points on 
v = Vrng ( r ) prior to p must lie to the future of points on 
V = vecif) prior to p, in the sense vrng i r ) > vedr) for 
r € (0,r*), where r* corresponds to p. Thus, the radial 
null geodesies must extend back to r = with v = vq 
satisfying v s < vo < Vah, and so must emerge from the 
singularity. 

It is noted that vah = w s = i>o = holds in the even- 
dimensional case. 



and 



dv 2 



(A22) 



in the untrapped and trapped regions, respectively 
Therefore, 



dr dr 



(A23) 



APPENDIX B: CURVATURE TENSORS 



We present the scalar curvatures in order to analyze 
the singularities. The curvature invariants for the metric 
(|3~7|) are 



if") 2 + ^^(/') 2 + 2) }° 3) (fc (Bl) 

o 2 . 

(B2) 



x 2 s Ra ^ = krr + D -±rr + ^(/o 2 - 2(J ' 2) f - 3 V - /)/' + (g - 2)( r 3) V - if. 

I r 2r r r 



I 

The Ricci scalar becomes These formulae are applicable not only for the DC-Vaidya 

solution but also for the DC-BTZ solution i|3.2fl or the 
_ 2) (D — 2)(D — 3) Vaidya-AdS solution (|C1|) . These quantities remain fi- 

R = —f / H 2 — f)' nite, as expected, for the three-dimensional BTZ solu- 



(B3) 



tion YTt 
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APPENDIX C: GENERAL RELATIVISTIC 
COLLAPSE 

In general relativity, the function f(v,r) in Eq. (|3.7|1 
and the energy density p for the D-dimcnsional Vaidya- 
AdS solution for k — 1 are 



/(„,r)-l-^ + rV 



and 



D-2 

2n 2 D r D - 2 



M, 



(CI) 



(C2) 



respectively, where M(v) is an arbitrary function of v. 
n 2 D is related to the /^-dimensional gravitational constant 
Go as k 2 d = 8itGd- We assume the form of M(v) as a 
power law 



M(v) = M v q 



(C3) 



as that in the case of DC gravity, where Mq(> 0) and 
q{> 1) are constants. We can see the central singularity 
appears at r = for w > 0. The singular nature at 
v = r = will be clarified below. The trapped region is 



M(v) = M v q > r u - A + I 



~ 2 r D ~ x . 



(C4) 



We will consider the D = 3 and D > 4 cases separately 
in the following. 



When Mo satisfies Eq. (IC9II for g = D — 3, the spacctime 
represents the formation of the naked singularities; oth- 
erwise we need a more careful investigation to determine 
the final state of the gravitational collapse. The right- 
hand side of Eq. I|C9(I goes to zero as D — ► +oo; here the 
naked singularity formation becomes less feasible as the 
spacetime dimensions are higher, similar to the case of 
the even-dimensional DC-Vaidya solution. While, there 
exist no null geodesies with the asymptotic form l|C5(l for 
1 < q < D — 3. The geodesies l|C6|) and (IC7|) are singular 
null geodesies for D — 3 < q < D — 1. Thus, the solution 
represents the formation of a massless ingoing null naked 
singularity for D — 3 < q < D — 1 for any Mq(> 0) and 
q = D - 3 with M satisfying Eq. (|C9|) . 

The asymptotic form of the null geodesic equation near 
v = r = is given by 







ith 



K°- 2 M (D - 3) 



S5 



(CIO) 



(Cll) 



for q = D — 3, where K$ satisfies the relation l|C8|) . This 
belongs to Case 1 in Sec. For g > D — 3, The asymp- 
totic form of the null geodesic equation near v = r = 
is given by 



r\2 



-jY^ + ser^+V) 







(C12) 



1. D > 4 case 



with 



From Eq. (|C4|) . wc can find that only the point v = 
r = has the possibility of being naked. Thus, a massive 
timelike naked singularity is absent in this case for any 
mass function with M(0) = and M > 0. 

We assume that the asymptotic form of the future- 
directed outgoing null geodesies near v = r = as 



v ~ K*r p , 



where K3 and p arc positive constants. After some 
straightforward calculations, we find the asymptotic so- 
lution 



v — 2r 



(C6) 



s 6 = 2"M q. 

This belongs to Case 2 with a = q — D 
It is then concluded 

A-0 * (1 + S 5 ) 2 



(C5) for q = D - 3, and 



hmA I) - (,J+1 V = S6(£'-2) 



(C13) 



(C14) 



(C15) 



for D — 3 < q < D — 1. The divergent behavior of tp is 
summarized in Table ITTT1 



for q > D — 3. For q = D — 3, which is valid for D > 3, 
we obtain 



t> ~ K 3 r, 



where K-i satisfies the relation 



(C7) 



TABLE III: The divergent behavior of ip around v = r = 
along the radial null geodesies with the asymptotic form (ICJ6I 
for D > q — 3 and jQZt f° r q = D — 3inD>4 dimensions. 



<2 


g = D-3 


D - 


-3<q<D- 


- 1 


q = D- 1 


9 > D-l 





A" 2 


■^-D + q + l 


const. 






M iff " 2 - A' 3 + 2 = 0. 



(C8) 



The condition for the existence for positive A3 satisfying 
Eq. (|C8|I is given by 



D — 3 case 



< M < 



1 



D-3 



D-2 \2(D - 2) 



D-3 



(Qg-\ From Eq. (|C4|I . the point r = 0, < v < vah = 

Mg -9 can be naked. In this case, the discussion in Sec. 



13 



IIVBI can be applied. There exist null geodesies which 
behave as Eq. (|4.14|l with n = 2 for < v < vah- The 
scalar curvature quantities are finite everywhere along 
these null geodesies although the energy density of the 
null dust fluid diverges. From Sec. IV Bl along the null 



geodesies from v = r = 0, ip = R^k^k 1 ' diverges as 
\~ 2+q . In contrast, along the null geodesies from r = 0, 
< v < vah, ip diverges as A -1 independent of the power 
of the mass function. 
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